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ABSTRACT 

The detailed abundance patterns of the stars within galaxies provide a unique window into the history of star 
formation (SF) at early times. Two widely used 'chronometers' include the a and iron-peak elements, which 
are created on short and long timescales, respectively. These two clocks have been widely used to estimate 
SF timescales from moderate-resolution spectra of early-type galaxies. Elements formed via s-process neutron 
captures (e.g., Sr and Ba) comprise a third type of chronometer, as the site of the main ^-process is believed 
to be intermediate and low-mass asymptotic giant branch stars. The [a/Ba] ratio in particular should provide 
a powerful new constraint on the SF histories of galaxies, in part because it is insensitive to the uncertain 
distribution of Type la SNe detonation times and the overall la rate. Here we present new measurements of 
the abundance of Sr and Ba in nearby early-type galaxies by applying stellar population synthesis tools to high 
S/N optical spectra. We find a strong anti-correlation between [Mg/Fe] and [Ba/Fe], and a strong positive 
correlation between [Mg/Ba] and galaxy velocity dispersion. These trends are consistent with the idea that 
more massive galaxies formed their stars on shorter timescales compared to less massive galaxies, and rule 
out several other proposed explanations for the observed super-solar [Mg/Fe] values in massive galaxies. In 
contrast, [Sr/Fe] ~ 0, with no strong variation across the sample. It is difficult to interpret the Sr trends without 
detailed chemical evolution models owing to the multiplicity of proposed nucleosynthetic sites for Sr. 

Subject headings: galaxies: stellar content — galaxies: abundances — galaxies: elliptical and lenticular, cD 



1. INTRODUCTION 

While it is clear that the majority of massive early-type 
galaxies formed the bulk of their stars at z > 1, the precise 
nature and duration of star formation in these galaxies has 
proven difficult to constrain. The difficulty arises from the 
fact that the main sequence turnoff point varies little with 
age at late times, implying that the integrated spectral energy 
distribution varies little with age at late times. For old stel- 
lar systems such as early-type galaxies, it has become com- 
mon to instead rely on elemental abundance ratios to probe 
the star formation histories at early times (e.g.. lWorthey et all 
Il992t Thomas et al. 2005). The classic example is the ratio 
of a to iron-peak elements, [a/Fe0 The a elements form 
mostly in massive stars, which evolve on short timescales, 
while iron-peak elements form mostly in Type la supernovae 
(SNe), which occur on longer timescales. The [a/Fe] ratio is 
thus a sensitive probe of SF on ~ 10 8-9 yr timescales (e.g., 
lTinslevlll979HThomas et al.l[l998h . 

One of the principal difficulties in using [a/Fe] to measure 
a SF timescale is the unknown Type la SNe delay time distri- 
bution. Other difficulties include the possibility that the over- 
all la rate varies in some systematic way with galaxy prop- 
erties, the possibility of selective mass-loss such that Fe is 
preferentially lost from th e system, and potent i al variation in 
the initial mass function (IWorthey et alJl!992t iThomas et all 
1999; Trager et al. 2000). A promising alternative chronome- 
ter is Ba, which is believed to form predominately within the 
envelopes of asymptotic giant branch stars via ^-process neu- 
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tron c aptures (Burbidg e et al.ll957l:lBu"sso et al.ll999l lHerwig 
l200l . The precise mass-dependent yields of Ba are presently 
not well-known, but there is hope that progress on this topic 
can be made in the near-term (in contrast to the delay time 
distribution for la's, which appears to be governed largely by 
the initial conditions of stellar binarity). The prospects for 
measuring Ba in the integrated light spectra of galaxies are 
favorable owing to the strong transitions of Ba II in the blue 
spectral region. 

Sr is another neutron capture element with strong transi- 
tions in the blue. Like Ba, Sr is predominantly produced 
by the ^-process, at least in the solar system. These two el- 
ements probe two of the three ^-process peaks, with Sr be- 
longing to the first (along with Y and Zr), and Ba belonging 
to the second (along with La, Ce, Pr, and Nd). While the 
nucleosynthetic origin of Ba is relatively secure, the same 
cannot be said for the elements in the first ^-process peak 
dCouch et alj|1974t IWooslev & Hoffman! 1 1992tlRaiteri et all 
l!993U Sneden et al. 2008). Indeed, current chemical evolution 
models of the Galaxy are unable to reproduce the observed be- 
havior of Sr-Y-Zr at low metallicity without appealing to ex- 
otic and/or ad hoc nucleosynth etic sites (e.g., Travagli o et alj 
12004 lOian & Wasserburg||2008l) . Perhaps the [Sr/Ba] ratios 
in metal-rich massive galaxies will provide some insight into 
this problem. 

In this Letter we employ our new stellar population syn- 
thesis (SPS) model to interpret high quality optical spectra of 
early-type galaxies in order to constrain the abundances of Sr 
and Ba. In Sections [2] and [3] we describe the model and data, 
and in Section|4]we present our results. 

2. MODEL 

The SPS model u sed he rein was developed in 
IConrov & van Dokkuml d2012al CvD12). The model 
adopts standard SPS techniques, including constructing 
libraries of isochrones and stellar spectra. Empirical spectra 
form the core of the model. The empirical stars are of approx- 
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Figure 1. Change in the spectrum due to a 0.3 dex enhancement in the [Sr/Fe] 
and [Ba/Fe] abundance ratios. The model is an integrated light spectrum with 
an age of 1 3 Gyr, solar metallicity, and a Kroupa ( 200 1 ) IMF. The spectra 
have been smoothed with a velocity dispersion of a = 350 kms -1 , which is 
the same smoothing scale applied to the SDSS data. The dotted line repre- 
sents the sensitivity limit for a signal-to-noise ratio of S/N= 10 3 , which is the 
typical S/N for the SDSS stacks in this wavelength range. 

imately solar metallicity and have solar abundance patterns. 
We computed a large grid of model stellar atmospheres and 
spectra in order to construct response functions (i.e., the 
relative change in the spectrum of a star due to a change in the 
abundance of a single element). The model atmospheres and 
spectra were computed with the ATLAS mod el atmosphere 
and spectrum synthesis package (iKuruczll 197011 19931) . ported 
to Linux by ISbordone et al J (120041) . Specifically, we use 
the ATLAS 12 code, computing new atmospheres for each 
change in abundance pattern. The line list was provided by R. 
Kuruc^E including linelists for TiO and H2O, amongst other 
molecules. The spectral response functions were applied to 
the empirical stellar spectra in order to create models with 
arbitrary abundance patterns. The model allows for arbitrary 
variation in the initial mass function (IMF) and spans ages 
from 3-13.5 Gyr. See CvD12 for further details regarding 
the model. 

We follow IConroy & van Dokkuml (12012b!) in fitting the 
model to data. In its present form the model contains 27 free 
parameters, including the redshift and velocity dispersion, a 
two-part power-law IMF, two population ages, four nuisance 
parameters, and the abundances of C, N, Na, Mg, Si, Ca, Ti, 
V, Cr, Mn, Fe, Co, Ni, Sr, Y, and Ba, and 0,Ne,S are var- 
ied in lock-step. These parameters are fit to the data via a 
Markov Chain Monte Carlo fitting technique. The data and 
models are split into four wavelength intervals (defined in the 
following section) and, within each interval the spectra are 
normalized by a high-order polynomial (with degree n where 
n = (A 

max 

A min )/100A; see iConrov & van Dokkuml 120 12b. 
for details). We have masked the spectral regions surrounding 
the Ha, [N II], [S II], H/3, [O III], and [N I] emission lines. 

The effective resolution of the data varies from R w 400 - 
1200 owing to intrinsic Doppler broadening. Both because of 
this low resolution and the fact that in our model the synthetic 
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Figure 2. Spectral indices for the stacked SDSS early-type galaxy spectra 
and the LRIS spectra of 35 individual galaxies. Index units are equivalent 
widths in A. The LRIS data do not cover the Ba6142 spectral feature. The 
Mg b Lick index is sensitive to Mg, while the o ther features are sensitive to 
Ba and Sr and are defined in Serven et al. (2005). 

spectra are only used differentially, several uncertainties asso- 
ciated with abundance analysis at high resolution (including 
the adopted microturbulent velocity, treatment of convection, 
non-LTE effects, and hyperfine splitting of certain levels) are 
mitigated to some extent herein. Moreover, Berg emann et al.l 
(120121) has demonstrated that the strong Sr II line at 4077A 
suffers from negligible non-LTE corrections at [Fe/H]« 0.0. 
Of course when working at low resolution one is more depen- 
dent on the fidelity of the model than the case where equiv- 
alent widths of individual unsaturated lines are directly mea- 
sured from the data. 

As an example of what can be measured from low resolu- 
tion spectra, in Figure [T] we show the change in the model 
spectrum due to a 0.3 dex increase in the [Sr/Fe] and [Ba/Fe ] 
abundance ratios (holding [Fe/H] fixed to the solar value). 
Th e reference mo del is for an age of 13 Gyr, with a Galac- 
tic jKr oupa 2001) IMF, and for solar metallicity. The model 
was convolved with a velocity dispersion of 350 kms -1 . No- 
tice that even at this low resolution, at least two Sr II lines and 
two Ba II lines produce a change in the spectrum of > 0.3% 
per 0.3 dex change in abundance. As we will demonstrate in 
SectionHl the existence of multiple, well-separated Sr and Ba 
lines allows us to test the robustness of our results in a rela- 
tively model-independent fashion. 

3. DATA 

In this Letter we consider two samples of early-type galax- 
ies. The first sample of 34 early-ty pe galaxies (and the nu- 
clear b ulge of M31) was presented in Ivan Dokkum & Conrovl 
d2012l) . The sample w as selected from t he SAURON sur- 
vey of nearby galaxies dBacon et al. 2001). High S/N spec- 
tra of these objects were obtained with the Low Resolution 
Imaging Spectrometer (LRIS) on the Keck I telescope. The 
wavelength coverage spans the range 3500A-10,400A with 
a gap at 5600A-7100A. Four wavelength intervals were de- 
fined, 0.4/im-0.46/im, 0.46/i/m-0.55/ito, 0.80/i/«-0.89/i»j, 
and 0.96/im- 1.02/im for fitting models to data. The wave- 
length choices reflect the limitations of the model and data 
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wavelength coverage and regions of severe telluric correc- 
tions. The S/N at 5000A varies from w 150 A" 1 to « 450 A" 1 . 
Spectra were extracted within R < R e /8, with a weighting in- 
tended to mimic a circular aperture, where R e is the effective 
radius. Th ese data have already been com pared to the CvD12 
models in IConrov & van Dokkuml d2012bl see the Appendix 
of that paper for a comparison of the best-fit models and data 
for each of the 35 LRIS galaxies). In that earlier work, we 
did not consider variation in Sr and Ba, and only briefly com- 
mented on the abundances of the other elements considered 
in the fit. We found that the ages, [Fe/H] and [Mg/Fe] values 
derived with our full spectrum fitting technique agreed with 
conventional Lick index-based techniques. 

The second source of data is provided b y the Sloan Dig- 
ital Sky Survey (SDSS; [York et all 120001). Following the 
methodology outlined in iGraves et al.l (120091) . we have se- 
lected a sample of galaxies within a narrow redshift interval 
of 0.02-0.06, with no detected emission in Ha nor in [O II], 
and with concentrated, de Vaucouleurs-like light profiles. The 
SDSS obtains spectra with fibers that have 3" diameters. The 
typical S/N of SDSS spectra in our sample is modest, of order 
20 A" 1 , and so we have chosen to stack the spectra in seven 
bins of velocity dispersion with mean dispersions in each bin 
of 88, 112, 138, 167, 203,246, and 300 kms" 1 . For the small- 
est sigma bin the SDSS fiber samples the inner 0.8 R e , while 
for the most massive bin the fiber samples the inner 0AR e . 
We have verified that our results do not change if we select 
galaxies within each bin such that the fiber samples the same 
fraction of R e . Each spectrum was continuum-normalized 
and convolved to an effective dispersion of 350 kms -1 be- 
fore stacking, and each spectrum contributed equally to the 
stack. The resulting S/N of the stacked spectra at 5000A 
ranges from w 500 A" 1 to 1800 A" 1 . As with the LRIS 
data, four wavelength intervals were defined for the SDSS 
stacks, 0.4/im-0.48/_wn, 0.48/_wn-0.58/jm, 0.58/-tm-0.64/j/n, 
and 0.80/ira-0.88/im for fitting models to data. These stacked 
data will be the focus of future work aimed at measuring the 
detailed abundance patterns as a function of galaxy proper- 
ties. In this work we limit our attention to the neutron capture 
elements Sr and Ba. 

As a first look at the data, in Figure we show spectral 
indices for both the SDSS stacked spectra and the individ- 
ual LRIS galaxies. The Mgb index is sensitive to Mg, while 
the other indices are sensitive to Ba and Sr lines (with central 
wavelengths indicate d by the index nam e). The Br and Sr in- 
dices are defined in Serve n et al. | (|2005h . while the Mg b index 
is defined in Worthe y et al.l (119941 no effort was made to place 
the latter on the Lick index scale). For this figure, the LRIS 
data were convolved to a common dispersion of 350 kms -1 
in order to afford a direct comparison with the SDSS stacked 
spectra. The LRIS data do not cover the 6142A spectral re- 
gion. Overall the agreement between the LRIS and SDSS 
samples is encouraging, given the very different observational 
setups and reduction techniques. It is however difficult to in- 
terpret spectral indices because an index is in general sensi- 
tive not only to the feature of interest but also to features con- 
tributing to the sidebands (i.e., the index measures a feature 
strength relative to a pseudocontinuum). This is especially 
true when the feature of interest is weak compared to other 
features. Full spectrum fitting such as the kind performed in 
the present work does not suffer from this drawback. 

4. RESULTS 
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Figure 3. Comparison between models and data in two wavelength regions. 
Upper panels show the continuum normalized spectra and lower panels show 
residuals in percent. The data (black lines) are for the stacked SDSS early- 
type galaxy spectrum in the dispersion bin centered on <r = 250 kms -1 . Two 
models are compared, one in which the [Ba/Fe] ratio is fixed to the solar value 
(red line), and a model in which the [Ba/Fe] value is included in the fit (with 
a resulting best-fit value of —0.4; blue line). The grey band denotes the S/N, 
and the dotted lines mark the strongest Ba II lines in this wavelength range. 
The thin line is the [Ba/Fe] response function (cf. Figure[TJ. Any single Ba II 
feature on its own is not particularly compelling, given the magnitude of the 
residuals elsewhere in the spectrum, but the coincidence of an improvement 
in x 2 around all four Ba II lines is strong evidence in favor of a sub-solar 
[Ba/Fe] abundance ratio in this particular spectrum. Moreover, the difference 
in Xmm between the two models is 1400, signaling strong preference for the 
modeiwith [Ba/Fe]= -0.4. 

Before showing the derived abundance trends, in Figure [3] 
we explore the quality of the fit around the Ba II lines for the 
stacked SDSS spectrum in the second highest a bin. In the fig- 
ure we compare two model fits, one in which [Ba/Fe]= +0.0 
and another in which the Ba abundance is allowed to vary in 
the fit (with a resulting best-fit value of [Ba/Fe]= -0.4). We 
also mark the locations of the strongest Ba II lines in this spec- 
tral region. It is noteworthy that the residuals decrease in all 
spectral regions sensitive to Ba. This is an important test of 
the results to follow because it is also evident from this figure 
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Figure 4. Derived [Ba/Fe] and [Sr/Fe] abundance ratios as a function of the 
derived [Mg/Fe] abundance ratio. Results are shown both for stacked SDSS 
early-type galaxy spectra and for the LRIS spectra of 35 individual galaxies. 
For the SDSS spectra, the mean velocity dispersion increases monotonically 
with increasing [Mg/Fe] . 

that there are numerous places where the residuals exceed the 
formal S/N limits of the data (demarcated by the grey band). 
If there were only one strong Ba line then it would have been 
difficult to argue that a region of high residual required non- 
solar [Ba/Fe] values. We return to this point later in this sec- 
tion. 

Our main result is shown in Figure^] In this figure we show 
the best-fit [Ba/Fe] and [Sr/Fe] abundance ratios as a function 
of [Mg/Fe]. We include both the stacked SDSS spectra and 
the LRIS spectra of individual galaxies. The formal statistical 
errors from the SDSS stacks are very small owing to the very 
high S/N of the spectra. It is encouraging that the LRIS data 
and SDSS stacks follow the same general trends of decreasing 
[Ba/Fe] and constant [Sr/Fe] with increasing [Mg/Fe]. Such 
low values of [Ba/Fe] are not uncommon amo ngst metal-poor 
halo stars in the Galaxy dSneden et al.ll2008l) . Indeed, below 
[Fe/H]« -3 essentially all halo stars have [Ba/Fe] < 0.0, with 
many stars having B a abundances as low as [Ba/Fe]«-2.0 
dFrancois et al.1 120071) . Such stars are also a enhanced, and 
have [Sr/Ba]> 0.0, again broadly consistent with the mean 
abundance ratios of the massive early-type galaxies (although 
with very different overall metallicities). 

Taken at face value, the [Ba/Fe] versus [Mg/Fe] trends are 
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Figure 5. Exploration of systematic effects on the [Sr/Fe] and [Ba/Fe] abun- 
dance trends. The fiducial model is compared to models in which only the 
blue (A < 5800A) spectral region is fit, in which there is no smoothing of the 
stacked spectra to a = 350 km s , and in which the IMF is fixed to the Galac- 
tic form. Results are shown only for the stacked SDSS early-type galaxy 
spectra for clarity. While these model variations induce changes of the order 
of 0.1 dex, the trends are clearly robust to these model permutations. 

consistent with the standard interpretation of [Mg/Fe] probing 
the SF timescale in these galaxies. A shorter timescale should 
lead to higher [Mg/Fe] and lower [Ba/Fe] values, as observed. 
Although not shown, the SDSS data also display a strong cor- 
relation between [Mg/Ba] and a. The derived [Fe/H] val- 
ues range from -0.1 to 0.0 over the SDSS early-type sam- 
ple, and so the trends observed in Figure |4] are unlikely to be 
due to metallicity-dependent yields. Other explanations for 
the [Mg/Fe]-er relation have been proposed, including a vary- 
ing IM F, varying SN la rate and selective mass-loss of SN la 
ejecta (IWorthey et al.l[l99l iThomas et al.lfl999t iTrager et all 
2000). The [Mg/Ba] data allow us to rule out the explanations 
that are particular to the la's, since the la's influence neither 
Mg nor Ba, but it cannot discriminate between a varying SF 
timescale and a varying high-mass IMF for the observed vari- 
ation in [Mg/Fe] with a. 

The near constancy of [Sr/Fe] is puzzling. It has been 
suggested that Sr may have multiple nucleosynthetic ori- 
gins, including the r-process, both the weak and main s- 
processes, and perhaps even hypernovae (iRaiteri et aL~ll!993b 
iTravaglio et"aTl 2004: Oi an & Wasserburgll2008l) . As fMg/Fel 
(and presumably the SF timescale) varies, then perhaps the 
main nucleosynthetic site varies in such a way as to keep a 
rough constancy in [Sr/Fe]. With detailed chemical evolu- 
tion models constrained to match the observed [Mg/Fe] and 
[Ba/Fe] trends, it may be possible to place interesting con- 
straints on the origin(s) of Sr in these massive galaxies. 

We turn now to several tests that have been performed to 
assess the robustness of these results. In Figure |5] we show 
the same abundance trends as before, for only the SDSS 
stacked spectra, and for a variety of permutations to the fidu- 
cial model. We cons ider a model i n which the IMF is held 
fixed to the Galactic (Kroupa 2001) form, a model in which 
only the blue (A < 5800A) spectral region is fit, and a model 
in which the SDSS stacks are not smoothed to a common res- 
olution of a = 350 kms~ . The latter is a particularly impor- 
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tant test because much more information is available in the 
low dispersion bins, which effectively have higher spectral 
resolution, compared to the high dispersion bins. Relatedly, 
when the data are not broadened to <j = 350 kms -1 , the asso- 
ciated response functions (Figure Q3 become much stronger. 
The derived abundance ratios vary between these model per- 
mutations by ~ 0.1 dex, but the overall qualitative trends are 
clearly robust to these details. We have also masked each of 
the two strongest Ba II lines (one at a time) and refit the mod- 
els to the data. The resulting [Ba/Fe] versus [Mg/Fe] trends 
are unaltered, although in the case of masking the Ba II 4555A 
line the overall [Ba/Fe] abundances shift lower by « 0.2 dex. 
This highlights that while trends appear to be robust in our 
analysis, there are zero-point uncertainties at the 0.1-0.2 dex 
level. 

A concern is that the features we are trying to measure are 
much weaker than other metal lines. In fact, as shown in Fig- 
ure |3]the expected Ba line strengths are similar to the typical 
systematic residual in the spectra after subtraction of our best- 
fit model! A robust test of our result is to randomly shuffle the 
central wavelengths of the strongest Ba II lines and recompute 
X 2 around the shuffled lines. In practice, this is achieved by 
altering the response function for Ba (see Figure 1) by assign- 
ing random central wavelengths in the interval 4000A-6000A 
while keeping the relative strengths and widths of the lines 
the same. We have performed this test 10 3 times for each of 
the SDSS stacked spectra. As an example of the results for 
the highest a bin, the residuals decrease around the strongest 
shuffled Ba II line in only 0.1% of the shuffled models when 
using the best-fit [Ba/Fe] value shown in Figure [4] More- 
over, in none of the shuffled models do the residuals decrease 
in both of the two strongest shuffled Ba II lines. In the four 
highest a bins, where the best-fit [Ba/Fe] ratios are < -0.1, 
the shuffled models are never able to reduce the residuals in 
all three of the strongest Ba II lines. It is therefore highly un- 
likely that the derived [Ba/Fe] abundances are the result of a 
chance alignment between the strong Ba II lines and spectral 
regions where the residuals happen to be high. 

The elemental abundances derived herein should provide 
novel constraints on the duration of star formation in early - 
type galaxies once chemical evolution models that include 
yields from the s- and r-processes are employed. The [Mg/Ba] 
ratio is particularly promising because the derived timescale 
will not depend on the uncertain delay time distribution or 
overall rate of Type la SNe, in contrast to timescales derived 
from the [Mg/Fe] ratio. In future work, we will present de- 
tailed abundance patterns for m 17 elements, which promise 
to provide further insight into the early formation histories of 
these now dormant galaxies. 
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